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Lepton-flavor-violating decays of the type L — > Z77 (fj, — » 677, r — > £77, and r — » ^77) are pro- 
posed as new probes of R-parity- violating supersymmetry. Non-penguin diagrams with a sneutrino 
C*") , that decays into two photons via a triangle graph might trigger such decays even in the absence of 

the corresponding radiative decays into one photon only, e.g. fi — > &y. Thus, processes of the type 
' L -— > £77 may provide an independent probe of new flavor physics. 

(N , 

[ PACS numbers: 12.60.Jv, 13.35.-r, 13.35.Bv 

Supersymmetry (SUSY) is one of the most promising candidates for new physics, curing some of the shortcomings 
of the Standard Model (SM), for instance, the emergence of quadratic divergences in the Higgs sector Q]. Despite its 
appealing theoretical features, however, there is so far no experimental evidence for SUSY up to the electroweak scale. 
Apart from direct evidence, i.e. direct production and decays of SUSY particles, indirect probes could be employed to 
search for SUSY, i.e. processes in which the SUSY particles emerge virtually either in loops or as tree-level mediators. 
\Q \ In this respect Flavor-Changing Neutral Current (FCNC) or Lepton-Flavor- Violating (LFV) processes are a natural 
searching ground for SUSY, since the SM background to such processes is either loop-suppressed (FCNC) or does not 
exist (LFV). 

In this paper the effects of new physics on a class of LFV processes, L — > Z77, are investigated, in which a lepton 
decays radiatively into a lighter one plus two photons. As an example, R-Parity- Violating (RPV) SUSY is discussed. 
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In the conventional RPV framework, the SUSY superpotential is supplemented with new RPV interaction terms 



luperp 

that are either baryon- or lepton- number- violating |2j, |3j. For the purpose of this paper, only the following lepton- 
i- £h ' number- violating operators are relevant: 

P 



S" 1 ! YJrpv 3 ^iiktabLIL^El + X' i . k e ab LfQ b j D'j c , (1) 



where Q and L are SU(2) doublet quark and lepton supermultiplets, respectively, and D c and E c denote the SU(2) 
singlet down-type quark and lepton supermultiplct; i, j,k — 1,2,3 label the generations and a, b are the SU(2) indices 
that are responsible for the antisymmetric property Ajjfc = —Xjik- 

Within the framework of RPV SUSY, /i — > e*y (and r — > ej, t — > fij) can be mediated by penguin-like diagrams of 
the type depicted in Fig. ^ Depending on the RPV couplings involved, the particles running in the loop are either 
quark-squark, or lepton-slepton. The process L — > Z7, where L and / are leptons of generations i and j, respectively, 
is then proportional to either \ imn Xjmn or A mni A mn j or A- m „A^- mn 0- 

So far, very little attention has been devoted to L — > Z77, under the assumption that it is always suppressed with 
respect to L — > Z7. An exception to this rule is discussed in |]|, where a leptoquark model in which L — > Z77 is much 
larger than L — > Ij is presented. Note that their model yields BR(/i — > 677) < 10 -18 ; see also 0. In general, an 
additional photon can easily be accommodated - one just has to add a photon to the diagram in Fig. ^ and obtain 
the X-reducible diagrams shown in Fig. (they are reducible in the sense that, after removing a photon line from 
L — > Z77, a legitimate RPV diagram for L — > Ij is obtained). Obviously processes mediated through such diagrams do 
not provide any further information on the corresponding RPV couplings, since they differ only by kinematic factors 
and by an extra electromagnetic coupling. 



*Work supported in part by the EC 5th framework, contract number HPMF-CT-2002-01663 
tphcga@physics.technion.ac.il 

shaouly@physics.technion.ac.il 
§ eilam@physics.technion.ac.il 
^Frank.Krauss@cern.ch 



FIG. 1: A penguin diagram responsible for the decays L 
is charged. RPV couplings are symbolized by a blob. 



?7 in RPV SUSY. The photon can couple to any line provided it 
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FIG. 2: The topology of (a) X-reducible and (b) X-irreducible diagrams for L — > £77. 



However, it is quite possible that the decay L — > l^ does not exist, or its rate is too small to be observed, whereas 
the decay L — > Z77 still has an appreciable rate. This could be the case if, for instance, all products of the type AA 
and A' A' are much smaller than products of the type AA', in which case the X-irreducible topology of Fig.JSJa may give 
rise to a much larger rate than the X-reducible diagrams of Fig. 

For example, if A122 and A 233 are the only non-vanishing lambdas, then — > ery) — while — * 677) =^ 0. 
Accordingly, such RPV coupling product combinations can be detected or constrained by the measured upper limit 
on the branching ratio of fx — > e77, which is at present 7.2 x 10 -11 @,@]. No experimental limits exist so far for the 
analogous decays of the r-lepton. 

There are proposals to measure (or constrain ) th e branching ratios of such SM-forbidden decays of the /i down to 
branching ratios of the order of 10~ 14 9], 10~ 15 |Toj or even down to 10 -18 This highlights the phenomenological 
impact of LFV muon decays. 

A general basis of effective operators describing the decays L — > ^77 with the topology of Fig. [3 i.e. L — > ISk — > I77, 
will be constructed first. Here Sk is a scalar particle of type k and i,j are generation indices of L,l, respectively. The 
momentum notation L(p) — > l(p')"f(k)"f(k') will be used. 

The amplitude can be defined as 



where m is the family index of the fermion that runs in the loop, and the index k corresponds to the intermediate 
scalar. Also, C^ ( - j7cm) , C^ ijkm) are Wilson-like coefficients with L(R) = (1 - (+)7 5 )/2, and 
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0\ = [ Ul (p')Lu L (p)} (k-k'g Xv -k' x k v )e* x e\ 
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O h 2 = [utip'^ULip)] ie^kXelel* (3) 

are effective operators. If the incoming lepton is right-handed, then L — > R in 0. 

In the framework of RPV SUSY with a sneutrino {pf.) exchange, the Wilson coefficients are given by: 

ol ^A* Afc 

C l#(ijkm)i C *2(ijkm)(t™ in thc loo P) = °! 4^ ' [f^M x )' 9lM%)] » 



j A* -A' 

C^fe-fcm); Ci*2(ijfem)( d m in the loop) = 0; N c Q 2 d - kjl fc " m ■ [f 1/2 (x); g 1/2 (x)] , (4) 



where Mz k is the mass of the A:th generation sneutrino and to/ , / = I or d, is the mass of the fermion of generation 
to that runs in the loop. Furthermore N c = 3, Qd = 1/3 and the abbreviation x = 2rnj/k ■ k' is being used. For 

diagrams with a charge- conjugate sneutrino (£/*) exchange, in (gj C^fefcm) <-> C^ 2 (jifem) ( i ' e - witn ^fcji ~^ ^feij) and 
M;j fc — > Mp*. The functions /]y 2 and g Y j2 (originating from the loop integral of the triangle) are given by 



1 + (1 — x) arcsin 2 
1 



1 
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f 1/2 (x) = -2x 

9\/2{x) — 2x arcsin 2 I —= I . (5) 



The emergence of these functions reflects the different structure of the integral related to the parity of the couplings 
of the fermion in the loop. The subscripts refer to the spin of the particle running in the triangle. Note that when 
x < 1, the arcsin (1/y/x) develops an imaginary part, which corresponds to the loop particle residing on its mass shell. 

In what follows, results will be given for both constituent and current light quark masses: rrid — 300 MeV, 
7ji s = 450 MeV, TO fc = 4.5 GeV, and m d = 10 MeV, m s = 120 MeV, m b = 4.2 GeV, respectively. We believe that 
using constituent quark masses makes more sense, since for our case here, it would allow the decay /1 — * edd, for 
current d quarks |12|. 

For the case of the muon decay, /1 — > e77, and for either down quarks with constituent masses or for the T-lepton 
running in the loop, the relation to 2 > to 2 > k ■ k' holds. Then, the limit to 2 ^> k ■ k' (to/ = m,d, to s , to^ or m T ) 
provides a good approximation. Therefore, thc loop integral functions can be replaced with their x — > 00 limits, 
fi/2(x ~ + 00) = — I and gi/2(x — * 00) = 2. Using this approximation, decay widths are given by 



,j f oi\ 2 1 Y&m' 
T(^e 77 )\ m}>k . k ,cA' RPV (-) — (6) 

where / = di (i = 1,2,3 for d, s or b, respectively) or / = r is thc fermion in the loop and S — Vk or z>£ {k = 1, 2 or 
3), is the scalar responsible for the flavor changing transition fj, — * e. Also 

A r R % =| A fcl2 A fc33 | 2 , A T R % =\ X k21 X k33 | 2 , A d R % =| X kl2 X' ku | 2 /9 , A d R % =| A fe21 A' te | 2 /9 . (7) 

The branching ratios scaled by the relevant RPV couplings for single particles running in the triangle are given in 
Table [I] In the spirit of this paper, only such down quarks and leptons running in the triangle are considered that 
would not give rise to the decay L — * Z7, i.e., only those leptons that cannot generate thc X- reducible-type diagrams. 

Comparing the current limit BR(fi — * 677) < 7.2 x 10~ n with the estimate for the maximal possible branching 
ratio for this decay mode in Table in conjunction with the existing bounds on the relevant RPV coupling products 
as given in Table ITT1 it can be argued that, at present, [i — > ejj does not impose any new constraints on the RPV 
SUSY parameter space. However, the situation will improve with the much more stringent experimental constraints 
that can be anticipated: on the third row of Table ITT1 the expected sensitivity of \i — > e77 to the relevant RPV coupling 
products is exhibited, assuming that future experiments will be sensitive to BR(/j, — > 677) 10~ 14 , see e.g. 0. 

For the case of r-decays, taking the existing bounds on the relevant RPV coupling products into account, the largest 
branching ratio is obtained for t — > ^77 with an s-quark in the loop. In particular, for | A233 A 32 2 1 = 0-01 (its upper 
bound, see [13), m s = 450 MeV and M 9 = 100 GeV, BR{t -> ^77) - 10~ 10 . This result is by far smaller than thc 
typical present limits on r decay branching ratios, which are of the order of 10 -7 , see e.g. |l ll . Moreover, planned 
experiments |15) . which will produce about 5 x 10 s taus per year, will also be insensitive to the decays r — > Z77. 

A few remarks are in order to conclude the discussion of the branching ratios: 
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Loop 

a- 


BR(fi -» e 77 )/|AAT 


BR{t -> e 77 )/|AAT 


BR{ T _> M77 )/|AAT 


(Afci2Afe mm ) 


(Afc2i Afc mm ) 


£fe 

(Afei3 A femm ) 


~* 

(Ami A fcmm ) 


(Afc23 Afc mm ) 


(Afc32 Afc mm ) 


d (m = 1) 


1.3 x 10" 8 (4.1 x 10" b ) 


8.7 x 10"' (5.6 x 10" 8 ) 


8.5 x 10"' (5.5 x 10" 8 ) 


s (m = 2) 


5.5 x 10" a (8.2 x 10~ 8 ) 


7.6 x 10"' (6.0 x 10"') 


7.4 x 10"' (5.9 x 10~ Y ) 


6 (m = 3) 


5.5 x 10" 11 (6.3 x 10" 11 ) 


2.8 x 10" 9 (3.2 x 10" y ) 


2.7 x 10" a (3.1 x 10" 9 ) 




Loop 
particle 

a- 


BR(n -> e 77 )/jAA|" 


BR(j -> e 77 )/|AA|" 


BR(t M77)/|AA|^ 


Vk 

(Afcl2 Afcmm) 


~* 

(Afc21 Afe mm ) 


(Afci3Afc mm ) 


~* 

(Afc3i Afe mm ) 


(Afc33Afe mm ) 


(Afe32 Afemm) 


e (m = 1) 


/i — > e 7 possible 


t — ► e 7 possible 


8.9 x 10" H 


/i (m = 2) 


/i — > e 7 possible 


5.0 x 10" b 


t — *■ ^i 7 possible 


r (m = 3) 


3.2 x 10" y 


t — > e 7 possible 


r — > /i 7 possible 



TABLE I: Branching ratios, scaled by the appropriate RPV couplings and for mp k = 100 GeV or = 100 GeV, for each loop 
particle. Results for quarks with both constituent and current (in parentheses) masses are given. Also indicated in this table 
are the relevant index combinations for the RPV coupling products when either a sneutrino or a charge-conjugate sneutrino is 
exchanged. 



RPV coupling 


Current limit 


Expected sensitivity from BR(fi — > e 77 ) < 10 14 


From I19| 


From 12, 20] 


A122A233 


none 


3 x 10" 3 c ' d 


1.8 x 10" 3 


A121 A133 


none 


2 x 10" 4 C ' E 


1.8 x 10" 3 


A122A211 


4.1 x 10" a 


4 x 10" 8 1 


8.8 x 10" 4 (4.9 x 10" b ) 


Al32A 311 


4.1 x 10" 9 


4 x 10" 8 1 


8.8 x 10" 4 (4.9 x 10" b ) 


A121 Am 


4.1 x 10" 9 


4 x 10" 8 1 


8.8 x 10" 4 (4.9 x 10" b ) 


A23lA 311 


4.1 x 10" 9 


4 x 10" 8 1 


8.8 x 10" 4 (4.9 x 10" b ) 


A122 A222 


7.7 x 10" 9 a >" 


9 x 10" 3 c ' 9 


1.4 x 10" 3 (3.5 x 10" 4 ) 


A132A322 


7.7 x 10" 9 a '° 


0.012 d ' h 


1.4 x 10" 3 (3.5 x 10~ 4 ) 


Al21 A 12 2 


7.7 x 10" 9 a ' b 


2.1 x 10~ 3 c 


1.4 x 10" 3 (3.5 x 10~ 4 ) 


A231 A322 


7.7 x 10" 9 a ' b 


0.012 d ' h 


1.4 x 10" 3 (3.5 x 10" 4 ) 


Al22 A2 3 3 


none 


7.4 x 10~ 3 c ' e 


0.014 (0.013) 


A132A333 


none 


0.016 d ' 1 


0.014 (0.013) 


Al2lA 133 


none 


6.9 x 10" b c ' e 


0.014 (0.013) 


A231A333 


none 


0.016 d '* 


0.014 (0.013) 



TABLE II: Current limits and expected sensitivity to RPV coupling products from future experiments. Limits are given for 
sparticle masses of 100 GeV. The limits shown on the fourth column scale as ms/100 [GeV], see text. 

"From (1 — » e conversion in nuclei. 

''This limit depends strongly on the strange-quark content in the nucleon and is, therefore, model-dependent. In particular, if the strange 
quark content in the nucleon is consistent with zero (see e.g. |2I| ). then this limit does not apply. 
c From charged current universality. 
d From T(r — » euu)/T(r — * {ivv). 
e From neutrino masses. 
^From {1 — » e conversion in nuclei |22| . 
"From BR(D+ -» K°* Vy)/ BR(D+ -* K°*e+v^). 
h From D° - D° mixing. Model-dependent. 
l Fromr w (Z°)/F ep (Z°). 



• There is also a second possibility for L — > Z 77 of the X-irreducible type to proceed: through the soft SUSY- 
breaking RPV terms [T(j CijkLiLjEf. or C^LiQjDl and slepton or squark in the loop. In this case the sneutrino 

couples to a pair of sleptons or down squarks of different helicity index, i.e. v/l/r, where / is a slepton or 
down squark, while the photons couple to fhfh or to Jr^r. Assuming for simplicity that the /l/r mixing is 
of O(l) and setting djk = MsusyCijk, C\- k = MsusyC 1 \jk, the width for L — ► Z 77 is (for rn 2 ~/m 2 L ^> 1 and 
mf/m 2 L « 1): 



V(T _» 7WI ~ 4« /S ( tt Y 1 Q ^ m L M SUSY m 

T ^^ l ^- A RPv{^) 256^8640 rrAm% ' (8) 
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where / = £ with Q i = —1 oi f = d with Qj = —1/3 if the sparticle in the loop is a slepton or down squark, 
respectively, and S = v or v* . Also, 



AtjtmVH —\r> \ i 2 A l]K - mU k —I C \ I 2 

RPV — I ^kmm^kij \ i ^RPV — I ^kmm^kji \ i 

lijd m u k _ q | f,, , i2 Aijdmv£ _ q I /}/ \ |2 [q\ 

RPV — y I ^ kram^kij \ i ^RPV — 3 | u kmra^kji \ ■ V y / 

Here m, fc = 1,2,3 denote the family of the slepton (£) or down squark (d) and sneutrino or antisneutrino, 
respectively. 

Plugging in numbers (e.g., setting M S usy = mg m = m 0k = 100 GeV), BR(^i -> 677) ~| C kmm \ki2 | 2 3.7x 10~ 14 
and BR(t — > 677) ~| Cfe m mAfci3 | 2 1.9 x 10 -12 . Both values are much smaller than the ones obtained with 
fcrmions in the loop. 

• If an RPV bilinear term exists, i.e. BLH U 16], then L — > Z77 can also proceed via the X-irreducible topology 
with one trilinear (A) and one bilinear (B) RPV insertion. Therefore, the decay would be oc A x B. In this 
case the off-shell sneutrino emitted from the decaying lepton mixes with the down-type Higgs boson, which then 
goes into the two photons with the well known (see e.g. 01) one-loop SUSY amplitude for H — ► 77 (with the 
replacement m 2 H — > 2k ■ fc'). This possibility will be investigated in detail elsewhere [l8|. 

• Another set of related processes that may be generated through the same X-irreducible diagrams are r — > (e, n)gg. 
Work on these processes is in progress [l8j. 

• The r-decays r — > (/1, 6)77 of the X-irreducible type that are considered in this paper, with the lepton loop, 
have tree-level "analogs" that probe identical combinations of lambdas: r — ► fiee and r — > /i^te [23| . However, 
the contribution from lepton loops is not necessarily the largest one: it was found above that the dominant 
contribution stems from the s-quark loop, for which the branching ratio can amount to up to O(10~ 10 ) (for 
r — > 1U77), taking into account the current constraints |13| . which are obtained from processes of the type 
r — > I + hadrons. 

• As it can be seen from the Table, fi—e conversion in nuclei is a strong rival to our processes. However, couplings 
that involve b quark and r lepton (and perhaps s quark) are beyond the reach of \x — e conversion experiments. 
Although our bounds on such processes are weak, they will strengthen if PRISM experiment will come into 
action. 

In this paper a new class of lepton flavor-violating-processes, L — > Z77 has been presented and discussed that might 
appear even if the corresponding decay into one photon only is suppressed or not existing at all. The underlying 
mechanism for such decays is the decay of the heavy lepton into the lighter one and a virtual scalar particle that 
goes into the two photons via a triangle diagram. Hence, such processes might probe the flavor structure of leptons 
coupling to these scalar particles that are a common feature of many models for new physics. Branching ratios for such 
decays have been calculated within the framework of R-parity-violating supersymmetry with sneutrinos as the scalar 
particles. Imposing RPV bounds emerging from other processes, the resulting branching ratios are well below current 
observation thresholds. However, for a planned new round of experiments, especially for fj, decays, these processes - if 
they exist - might shed light on potential new physics or in turn might help to set more stringent exclusion bounds. 
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